APS/123-QED 



o 

(N 



Coherent and incoherent processes responsible for various characteristics of nonlinear 

magneto-optical signals in rubidium atoms 
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We present the results of an investigation of the different physical processes that influence the 
shape of the nonlinear magneto-optical signals both at small magnetic field values (~ 100 mG) and 
at large magnetic field values (several tens of Gauss). We used a theoretical model that provided 
an accurate description of experimental signals for a wide range of experimental parameters. By 
turning various effects "on" or "off" inside this model, we investigated the origin of different features 
of the measured signals. We confirmed that the narrowest structures, with widths on the order of 
100 mG, are related mostly to coherences among ground-state magnetic sublevels. The shape of 
the curves at other scales could be explained by taking into account the different velocity groups 
of atoms that come into and out of resonance with the exciting laser field. Coherent effects in the 
excited state can also play a role, although they mostly affect the polarization components of the 
fluorescence. The results of theoretical calculations are compared with experimental measurements 
of laser induced fluorescence from the D2 line of atomic rubidium as a function of magnetic field. 
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PACS numbers: 



I. INTRODUCTION 



When coherent radiation excites an atomic system 
with ground-state angular momentum F g and excited- 
state angular momentum F ei coherences can be created 
among the magnetic sublevels [HE]. At low laser in- 
tensity, coherences appear in the excited state of the 
atom. As the laser intensity increases, the absorption 
processes become nonlinear, and coherences are created 
among the magnetic sublevels of the ground state as 
well. When the degeneracy among the magnetic sub- 
levels is lifted by applying an external field (in our case 
magnetic), the coherences are destroyed. As a result, 
nonlinear magneto-optical resonances can be observed in 
the laser-induced fluorescence (LIF) plotted as a function 
of magnetic field. For linearly polarized radiation excit- 



ing a transition F g 



F R = F a 



1, these resonances 



will be bright, that is, the atoms will be more absorbing 
at zero magnetic field (3HI>]. When F e < F g , the reso- 
nances will be dark, or less absorbing at zero magnetic 
field [8]. The NMOR features can be as narrow as 
1CU 6 — 1CU 5 G when buffer gas or antirelaxation coating 
of the cell is used because of the slow relaxation rate of 
the ground state [5] . This characteristic makes them suit- 
able for many applications, such as, for example, mag- 
netometry [TO], lasing without inversion [TT], electrically 
induced transparency [T2] , slow light and optical informa- 
tion storage [131 E] , atomic clocks [TS] , and narrow-band 
optical filters [16j . However, these narrow resonances are 
usually found within broader structures with features on 
the order of several Gauss or several tens of Gauss in a 
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plot of LIF versus magnetic field. Our study focuses on 
these broader structures, which are interesting in them- 
selves and also for some practical applications at higher 
magnetic field values, like optical isolators [17]. Using a 
theoretical model that has been developed over time and 
mostly was used to describe the narrow magneto-optical 
resonances but can reproduce the magneto-optical sig- 
nals with high accuracy over a large range of magnetic 
field values |18| . we investigated the peculiar shape and 
sign (bright or dark) of these structures, as well as the 
physical processes that give rise to them. 

In order to describe magneto-optical signals over a 
magnetic field range of several tens of Gauss or more, 
it is necessary to include in the model excited-state co- 
herences, energy shifts of the magnetic sublevels in ex- 
ternal fields, which bring levels out of resonance with 
the narrow-linewidth laser radiation, and the magnetic- 
field-induced mixing of the atomic wavefunctions, which 
changes the transition probabilities of the different transi- 
tions between ground and excited-state sublevels [T9"ll20| . 
Moreover, it is necessary to treat various relaxation pro- 
cesses, the coherence properties of the laser radiation, 
and the Doppler effect. Since at least the 1970s, magneto- 
optical signals in alkali atoms have been modeled by 
solving the optical Bloch equations for the density ma- 
trix [2T]. Simple models were able to describe the nar- 
row resonances fairly well [22], but failed to describe the 
signals at fields of several Gauss or more. With time, 
these models become more sophisticated as the afore- 
mentioned effects were incorporated [23-25J, and now 
the agreement is often excellent, at least up to magnetic 
fields over one hundred Gauss. Thus, numerical models 
have become useful tools for understanding the physical 
processes that give rise to various features in the signals, 
because different physical processes can be included in 
the models or excluded one by one. Analytical studies, 



on the other hand, can demonstrate more explicitly a link 
between a particular physical process and the observable 
outcome. Thus, in [35] analytical formulae were devel- 
oped that allow one to calculate the contrast of bright 
resonances. In another study, a theoretical model of the 
electromagnetically induced absorption (EIA) was con- 
structed for a hypothetical F g = 1 — > F e = 2 transi- 
tion [35]. It was possible to show from a purely theo- 
retical point of view that the sub-natural linewidth res- 
onance in EIA was related to the transfer of coherence 
from the excited state to the ground state. More recently, 
sophisticated analytical models were developed that are 
valid in the low-power region, and were applied to ex- 
perimental measurements on the cesium Dx line j27[ [28] . 
Comparison with experiments confirmed that the narrow 
resonances arise when polarization is transferred from the 
excited state to the ground state. In [29] an analytical 
model was used to analyze the influence of partially re- 
solved hyperfine structure in the ground or excited state 
on nonlinear magneto-optical rotation signals. Numerical 
studies such as ours can complement these analytical in- 
vestigations, because the numerical models can be made 
to apply over a wider range of laser power densities and 
consider realistic, Doppler-broadened atomic transitions 
in the manifold of the hyperfine levels, that is to say, take 
into account multiple adjacent transitions. 

Our study focused on the D 2 line of 87 Rb as a model 
system. Since the origin of the narrow structure had al- 
ready been shown to be connected to coherences in the 
ground state [23I2E]: our study primarily aimed at under- 
standing the wider features of the magneto-optical signals 
up to magnetic field values of several tens of Gauss, since 
such understanding is important in itself and will help 
to improve the models of the narrow resonances used for 
applications. Nevertheless, since a numerical model such 
as ours gives complete flexibility to turn different effects 
"on" and "off", we were also able to confirm the origin 
of the narrow structure using a different technique, i. e., 
one that is not analytical. 

The level structure of the transition studied here is 
shown in Fig. n] 30] . The transition was excited by lin- 
early polarized laser radiation. Figure [2] shows the rel- 
ative transition probabilities from the ground-state sub- 
levels of the F g = 2 level to the excited-state sublevels of 
the F e = 3 level when the linearly polarized exciting ra- 
diation is decomposed into coherent circularly polarized 
components. It is assumed that the light is polarized 
perpendicularly to the direction of the external magnetic 
field (see Fig. [3J) This scheme implies that Am = 2 co- 
herences are created between different Zeeman sublevels 
in the excited state as well as in the ground state. Two 
distinct processes contribute to ground-state coherence. 
The first process creates coherence in the ground state 
through direct interaction with the radiation field via 
A-type absorption. In the second proccess the V-type 
absorption creates coherences in the excited state, which 
then can be transferred back to the ground state via spon- 
taneous emission, see Eq. (13.13) in [2J. Fig. [2] shows 
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FIG. 1: Scheme of the hyperfine levels and allowed transi- 
tions of the Z?2 line of 87 Rb. 
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FIG. 2: Relative transition strengths from the ground- 
state magnetic sublevels to the excited-state magnetic sub- 
levels when the linearly polarized exciting radiation is de- 
composed into a circularly polarized components for the 
F g — 2 — s> F e = 3 transition of the D2 line. The Lande 
factor gp is given at the left of each particular hyperfine level 



that both V-type and A-type transitions are present in 
our physical system. 

The paper is organized as follows: Sec. II outlines the 
theoretical model. In Sec. Ill we describe the experi- 
mental conditions, and in Sec. IV we discuss the results 
and attempt to decompose the modeled signal into com- 
ponents that are related to different physical processes. 



II. THEORETICAL MODEL 

The theoretical model is based on the density ma- 
trix approach. The density matrices are written in the 
\^ 1 Fi 1 mp) basis where Fi denotes the quantum number 
of the total atomic angular momentum, mp, the respec- 
tive magnetic quantum number, and £, all other quantum 
numbers. The time evolution of the density matrix is de- 
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FIG. 3: Geometry of the excitation and observation direc- 
tions. 



scribed by the optical Bloch equations [3T] 



ih 



dp 

at 



H,p 



ihRp, 



(1) 



which include the full atomic Hamiltonian H = Hq + 
Hb + V constructed from the unperturbed atom's Hamil- 
tonian Hq , which depends on the internal dynamics of the 
atom, the Hamiltonian Hb, which describes the atom's 
interaction with the external magnetic field, and the 
dipole operator V, which represents the atom's inter- 
action with the electromagnetic radiation. The inter- 
action with the magnetic field gradually decouples the 
total electronic angular momentum J and nuclear spin I, 
which means that F no longer is a good quantum num- 
ber, while m still remains a good quantum number. To 
deal with this effect, mixing coefficients between different 
hyperfine states in the magnetic field are introduced in 
the model. The relaxation operator R in (fiT) accounts for 
the spontaneous decay that transfers atoms from the ex- 
cited state to the ground state, the collisional relaxation, 
and the transit relaxation. The latter occurs when atoms 
leave and enter the interaction region as a result of their 
thermal motion. 

The optical Bloch equations can be written explicitly 
for each element of the density matrix. Applying the 
rotating wave approximation and assuming the density 
matrices do not follow promptly the random phase fluc- 
tuations of the electromagnetic radiation, we may decor- 
relate the time-dependent differential equations from the 
fluctuating phase and average over it. Thus we may adi- 
abatically eliminate the equations that describe the opti- 
cal coherences and obtain rate equations for the Zeeman 
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In both equations of Q the first term describes the opti- 
cally induced transitions to the level described by a par- 
ticular equation, and the second term, the transitions 
away from it, with dij being the element of the dipole 
transition matrix that can be calculated according to the 
Wigner-Eckart theorem [5] . The terms S g . e . and complex 
conjugate H*. ff are described below. The third term de- 
scribes the coherence destruction by the magnetic field, 
with wy = ' > J denoting the energy difference between 
levels \i) and \j) caused by both the hyperfine splitting 
and the nonlinear Zeeman effect. The fourth term de- 
scribes relaxation due to transit relaxation, collisions and 
spontaneous decay (only for the excited state). Two ad- 
ditional terms in ( |2a[ ) stand for population transfer to 
the ground state via spontaneous decay from the excited 
state (fifth term) and unpolarized atoms entering the in- 
teraction region as a result of their thermal motion (sixth 
term). 

The symbol S g . ej in equation ^ describes the strength 
of interaction between the laser radiation and the atoms 
and is expressed as follows: 



fl 2 



r+7+Ac 



+ i(u — k^v - 



(3) 



J gi ej 



where flu is the Rabi frequency, further discussed in 



Sec. [TV) T and 7 are the rates of spontaneous decay and 
transit relaxation, Aw is the finite spectral width of the 
exciting radiation, u> is the central frequency of the ex- 
citing radiation, k.Q the respective wave vector, and k^v 
is the Doppler shift experienced by an atom moving with 
a velocity v. The dependence of the absolute value of 
S g . e . at fixed i and j on the magnetic field is responsible 
for the effects of magnetic scanning discussed in Sec. |IV[ 
while the imaginary part of S 9ie . represents the dynamic 
Stark effect. 

The steady state solution of the rate equations (pi) 
yields the density matrices that describe population of 



magnetic sublevels and Zeeman coherences of both the 
ground and excited states. The density matrix of the ex- 
cited state is used to calculate the fluorescence signal for 
an arbitrary polarization component e: 



Ifi(e)=Io 



9i,ej,e k 



j*(ob) i{ob) 



ft e i 



efc3i 



(4) 



(ob) 



where Iq is a proportionality coefficient and d ejgi are 
elements of the dipole transition matrix for the chosen 
observation component e. The unpolarized fluorescence 
signal in a particular direction was calculated by sum- 
ming over two orthogonal polarization components. To 
take into account the Doppler effect, this quantity was 
averaged over the one-dimensional Maxwellian distribu- 
tion of atomic velocity along the direction of the laser 
beam propagation axis. In addition, the density matri- 
ces for some particular velocity groups are used to obtain 
angular momentum probability surfaces [H [3H [53] . 



III. EXPERIMENT 

The experiments were carried out at room temperature 
on natural mixture of rubidium isotopes in a cylindrical 
Pyrex vapor cell with optical quality windows, 25 mm 
long and 25 mm in diameter, produced by Toptica, A.G. 
of Graefelfing, Germany. The geometry of the excita- 
tion and observation is shown in Fig. [3j The 780 nm 
exciting laser radiation propagates along the x axis with 
linear polarization vector E pointing along the y axis. 
The total LIF (without polarization or frequency dis- 
crimination) was observed along the z axis, which was 
parallel to the magnetic field vector B. The laser was 
a homemade extended-cavity diode laser. The magnetic 
field was supplied by a Helmholtz coil and its value was 
scanned by controlling the current in a Kepco BOP-50-8- 
M bipolar power supply. Signals were recorded by a pho- 
todiode (Thorlabs FDS-100). The laser frequency was 
determined by means of a saturated spectroscopy setup 
in conjunction with a wavemeter (WS-7 made by High- 
Finnesse). The beam profile was measured by means of 
a beam profiler (Thorlabs BP104-VIS). The full width 
at half maximum was assumed to be the beam diame- 



IV . The 



ter used in the calculations [see Eq. (|6j), Sec. 
ambient magnetic field along the x and y directions was 
compensated by a pair of Helmholtz coils. The entire ex- 
perimental setup was located on a nonmagnetic optical 
table. 



IV. RESULTS AND DISCUSSION 

As the main tool for our present investigation was a 
numerical model, the first step was to show that it accu- 
rately described the measured signals over a large range 
of magnetic field values. Previous studies had already 



shown the model to be accurate in many experimental sit- 
uations |34] in which narrow magneto-optical resonances 
form in weak magnetic fields (B < 0.3 G) as a result 
of coherences created among the magnetic sublevels of 
the ground state. Figure [4] shows plots of LIF versus 
magnetic field over the range —40 G to +40 G when the 
laser was tuned to the F g = 2 — > F e — 3 transition of 
the £>2 line of 87 Rb at different laser power densities, as 
well as a plot of contrast versus laser power density. It 
must be noted that due to proximity of other hyperfine 
levels in the excited state, the Doppler effect and mag- 
netic scanning, other hyperfine levels were also excited 
at least partially. These transitions are included in our 
theoretical model as well. We defined the signal contrast 
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where I m in is the minimum LIF value (zero first deriva- 
tive and positive second derivative) around B = 0, and 
Imax is the LIF value at the first point with vanishing 
first derivative and \B\ > 1 G. Filled circles represent 
experimentally measured values, whereas the line shows 
the result of a theoretical calculation. In order to obtain 
an appropriate fit to the data, it was necessary to adjust 
two parameters. The first parameter was the constant 
fc 7 that relates the ratio of the mean thermal velocity Vth 
of the atoms and the characteristic diameter of the laser 
beam d to the transit relaxation rate 7 as 



7 = k 



Vth 

7 d 



Icol 



Vth 



(6) 



where 7 co ; is the rate of inelastic atom-atom collisions. 
An estimated value of j co i at room temperature, assum- 
ing the spin-exchange cross section for Rb-Rb collisions 
a sw 2 • 10~ 14 cm 2 [35] is several orders of magnitude less 
than the first term in Eq. (pi), and so it was omitted 
in the actual calculations. The second parameter Icr re- 
lated the Rabi frequency Ojj to the square root of the 
the experimental laser power density I according to 



Qfl = k T 



(7) 



where ||d|| is the reduced dipole matrix element that re- 
mains unchanged for all transitions within the D 2 line 
at a well documented value 2 , and c is the speed of 
light. Both fitting parameters (fc 7 and kji) would be 
equal to unity for a rectangular beam profile of the ex- 
citing laser and atoms moving with the mean thermal 
velocity across the middle of the beam profile. In our 
experiment the beam profile is roughly Gaussian, and so 
the laser beam diameter cannot be defined unambigu- 
ously. Furthermore, atoms are moving along random 
trajectories with velocities distributed according to the 
Maxwellian velocity distribution. Thus we allow the val- 
ues of these constants to deviate from unity in order to 
obtain an optimal fit between the modeled and experi- 
mentally recorded results. 



The actual values of the fitting parameters were fc 7 = 
0.5 and ku — O.fi. These values indicate that the 
interaction of atoms and laser radiation in the wings 
of the (roughly Gaussian) beam profile cannot be ne- 
glected. Thus for a beam with d — 1.6 mm (estimated 
in the experiment as defined in Sec. Ill I and laser power 
P = 20 /iW, we obtained the following values that were 
used in the modeling: 7 = 95 kHz and FIr — 0.75 MHz. 
Another important parameter for modeling and inter- 
preting the results is the natural linewidth, which is 
r = 6.067 MHz [3D]. Having obtained the optimum val- 
ues for these parameters by trial and error, these values 
were used to fit simultaneously all experimental data ob- 
tained for different transitions and different values of the 
laser power density. (The top left plot in Fig. |4]was mea- 
sured in a different experiment dedicated to the narrow 
structure [34] , and so the experimental conditions and fit- 
ting parameters were slightly different in this case, and 
the range of the measured magnetic field was smaller.) 
Agreement between experiment and theory was rather 
satisfactory, which shows that the model serves as a good 
basis for understanding the dependence of LIF on the 
magnetic field over a broad range of magnetic field val- 
ues. 

The narrow resonance at zero magnetic field is related 
to the destruction of coherences in the ground state by 
the magnetic field as we will show in the next paragraphs. 
Under our experimental conditions, it had a width of 
about one hundred milligauss and was clearly visible right 
at zero magnetic field. A detailed study of this reso- 
nance was performed in |34| showing that this structure 
points up or down (changes the sign of the second deriva- 
tive) depending of the laser power density. Under the 
present experimental conditions it appeared as a narrow 
structure with a negative second derivative (pointing up- 
wards) . This narrow resonance was located in the center 
of another structure with a positive second derivative and 
a width of several Gauss. 

In order to study how different physical effects influ- 
ence those features of the signal that appear at different 
scales of the magnetic field, we used the same theoreti- 
cal model, while turning different physical processes "on" 
and "off". Three processes were considered: destruction 
of ground-state coherences by the magnetic field, destruc- 
tion of excited-state coherences by the magnetic field, 
and the "Zeeman magnetic scanning effect" , which in- 
volved optical transitions between different Zeeman sub- 
levels that come into resonance with the laser radiation as 
a function of the magnetic field strength and the atomic 
velocity. The results are shown in Fig. [5] When all effects 
were included, we obtained structures on the scale of 100 
mG, several Gauss, and several tens of Gauss [see Fig. 
[51(a)]. The latter two are, as we will show later, caused 
by detuning effects as the hyperfine levels are split in the 
external magnetic field, and we will refer to these features 
as the "wide structure" . 

When the effect of the changing magnetic field on the 
coherences was neglected, which was done by setting the 



third term in Eq. pi) to zero for both ground and ex- 
cited states [Fig. [5p(b)], the small, narrow peaks dis- 
appeared completely, whereas the other structures re- 
mained largely, but not completely, unchanged. In order 
to consider only the ground-state coherence effects, the 
excited-state coherences were decoupled from the mag- 
netic field by setting the third term of Eq. ( 2b I to zero, 
and the detuning effects were "turned off" by taking the 
term 0J giej in the denominator of (pf to be independent 
of the magnetic field and keeping its value at the value 
it has at B = 0. Only the narrow structure was repro- 
duced when only the magnetic field's destruction of the 
ground-state coherence effects were taken into account 
in Fig. [5] (c) . The results shown in Figures [5] (b) and 
pKc) clearly attributes the narrow structure to the ground 
state coherences and their destruction by the magnetic 
field. At the same time the resonance with a width of 
several Gauss in Fig. [5] (b) is seen to be related to de- 
tuning effects, which where the only ones considered in 
that calculation. 

When only the excited-state coherent effects were 
taken into account in a similar way, a structure with 
negative second derivative and a width of several Gauss 
appeared; the contrast was only one or two percent [Fig. 
[5] (d)]. The structure had the same characteristic width 
(T « WAm=2) as the linear Hanle effect of the excited 
state [36]. The linear Hanle effect cannot be observed in 
our experiment as it requires discrimination of the po- 
larization components of the LIF, and so we attribute 
this structure to the nonlinear Hanle effect of the excited 
state. Calculations at several Rabi frequencies showed 
that the peak associated with this effect became smaller 
as the Rabi frequency changed from 1.0 MHz to 2.0 MHz. 
Moreover, at 2.0 MHz another small dip with positive sec- 
ond derivative appeared inside the peak at zero magnetic 
field; a further increase in Rabi frequency indicates a sim- 
ilar behaviour though on different scale as in Fig. [5] (c) 
(the effects produced by the destruction of ground state 
coherences). In any case, the calculations show that ex- 
cited state coherences play no role in the narrow struc- 
ture. 

The main origin of the wide structure can be under- 
stood by considering Fig. [6] The left panel shows how 
a magneto-optical signal can be decomposed into contri- 
butions from different velocity groups. The solid black 
line represents the signal of a vapor at room tempera- 
ture and is formed from an average over all the veloc- 
ity groups in the Doppler profile. The dashed and dot- 
ted lines represent contributions from different velocity 
groups. One can see that the superposition of the con- 
tributions from the dashed and dotted lines would yield 
a shape similar to the black line. The right panel ex- 
plains why each velocity group has its own shape. The 
laser is assumed to be on resonance at zero magnetic field 
with a group of atoms that is stationary with respect to 
the propagation direction of the laser radiation (v x = 0) 



for the F„ 



F e — 3 transition. All other velocity 



groups therefore interact negligibly with a laser field that 
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FIG. 4: (Color online) LIF versus magnetic field value for the F g — 2 

laser power density 7: (a) 0.14 mW/cm 2 , (b) 1 mW/cm 2 , (c) 10 mW/cm 2 . The bottom right panel shows the contrast of the 
central minimum as a function of laser power density. Filled circles correspond to experimentally measured values, whereas the 
solid line shows the result of a calculation. Note the different scales in (a) and (b-c). 



is detuned by the Doppler shift. As the magnetic field is 
applied, all magnetic sublevels shown in Fig. [2j except 
those with m = 0, are shifted as a result of the Zeeman ef- 
fect. We may say that a magnetic scanning is performed 
by bringing into resonance a group of atoms with some 
velocity v x — v(B). The function v(B) in general is non- 
linear and is explicitly determined by the nature of the 
(nonlinear) Zeeman effect. As a result of the magnetic 
scanning, the shapes of the angular momentum distribu- 
tions induced by the laser radiation differ as a function 
of magnetic field for each velocity group, which can be 
explicitly shown by the angular momentum probability 
surfaces [32 [33] for the excited state. When the angu- 
lar momentum probability surfaces are drawn, only the 
F e — 3 hyperfine level is taken into account, as other 
hyperfine levels are far away from resonance for the mag- 
netic field values and velocity groups shown in Fig. [6j and 



their populations are negligible. We may anticipate from 
Fig. [Gland the preceding discussion that, at a particular 
magnetic field value, some group of atoms with corre- 
sponding velocities becomes effectively oriented in either 
the positive or the negative direction of the axis along 
which the magnetic field is applied. Further, the whole 
ensemble of atoms becomes aligned along the same axis 
at magnetic values that produce the LIF maxima around 
±10 Gauss. 



V. CONCLUSION 

Nonlinear magneto-optical resonances from the Z?2 line 
of 87 Rb have been studied experimentally and theoreti- 
cally up to magnetic field values of 40 G. The theoretical 
model was based on the optical Bloch equations and in- 
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FIG. 5: (color online) Theoretical calculations of LIF versus magnetic field B for the F g — 1 — > F e — 3 transition of 87 Rb 
with different physical effects taken into account: (a) all effects taken into account, (b) detuning effects only, (c) ground state 
coherence effects only, (d) excited state coherent effects only. Note the different scales! The parameters used in the simulation 
were as follows: 7 = 0.019 MHz, Auj Laser = 2 MHz, a Dopp i er = 216 MHz, D S te P ~ 1.73 MHz 



eluded the coherence properties of the laser radiation, all 
adjacent hyperfine transitions, the mixing of magnetic 
sublevels in the external magnetic field, and the Doppler 
effect. The model described the experimentally measured 
signals very well. By removing individual physical pro- 
cesses from the model, it was possible to deduce the phys- 
ical origin of the different features observed in the signals. 
As expected, the narrow structure was related to coher- 
ences among ground-state Zeeman sublevels induced by 
the exciting laser radiation. Coherences among excited- 
state sublevels were found to have a small effect on signals 
at magnetic field scales of several Gauss. The origin of 
the wide structure was explained in terms of contribu- 
tions from different velocity groups. With these results, 
it is possible to understand the origin of the variation in 
LIF as a function of magnetic fields in the range up to at 
least several tens of Gauss. 



size the necessity to incorporate a number of processes 
in a theoretical model that aims to provide a quantita- 
tive description of magneto-optical effects. The most im- 
portant of these effects are 1) the Doppler effect, 2) the 
magnetic scanning and 3) the change in the transition 
probabilities due to the magnetic mixing of the hyper- 
fine levels, which can reach 30% for 87 Rb D2 excitation 
at B = 40 G. Although each of the processes can be 
treated separately to obtain an analytical description, in 
order to have an accurate description that is valid over 
a wider range of laser power densities and magnetic field 
values one has to treat all the processes simultaneously. 
On the other hand, a numerical model that incorporates 
a number of processes can be used to estimate limiting 
conditions for various approximations used in analytical 
models in the way described above. 



We may conclude that the results of this study empha- 
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FIG. 6: Decomposition of a magneto-optical signal into a superposition of signals from different velocity groups and at 
different magnetic fields. Left panel: The solid, black line shows the magneto-optical signal as it would be observed in a vapor 
cell at room temperature. The dashed and dotted lines show the signals for the different velocity groups that make up the 
room temperature velocity distribution. Right panel: Distribution of the atomic angular momentum at different values of the 
magnetic field B for the velocity groups in resonance at a (Doppler) detuning of MHz, 5 MHz, and -5MHz. 
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